日本沼虾微卫星标记开发、遗传结构分析及性别相关遗传标记筛选研究

【摘要】日本沼虾隶属十足目、长臂虾科、沼虾属，自然分布于中国、日本、越南及俄罗斯远东地区，在我国各地淡水水域中都有分布。日本沼虾的养殖始于二十世纪五十年代，从2004年起，我国日本沼虾的年养殖产量就超过了20万吨。多年来，日本沼虾在市场上一直供不应求，畅销不衰，是一种养殖效益非常可观的淡水虾类。但长期以来，繁殖日本沼虾用的亲本多为野生个体，未经过系统选育，日本沼虾的养殖单产较低，且多是套养在各种鱼类的养殖池塘中，主养日本沼虾的池塘较少。因此开展日本沼虾的种质资源保护、良种选育工作是当前日本沼虾养殖业中亟待解决的问题。由于当前已知的日本沼虾微卫星标记只有二十多个，这严重地限制了日本沼虾的群体遗传结构分析、种质资源评价等工作；此外，日本沼虾雄性生长速度显著大于雌性，开展单性苗种培育意义重大，但当前没有鉴别日本沼虾遗传性别的分子标记，使得单性苗种培育工作一直无法得以开展。本论文的研究内容主要包含如下3个方面：1.日本沼虾多态性微卫星标记的开发：（1）FIASCO法开发微卫星标记。通过生物素标记的探针（GT）13和（GA）13富集微卫星序列，共获得含微卫星的序列334个，成功设计引物110对，设计引物的序列含有的微卫星重复多为2碱基重复，重复数多数在10次以上。设计的引物扩增的目的条带大小在95-356bp之间。多态性分析显示，24对引物扩增出的条带为清晰稳定的多态，开发效率（多态性引物占总设计引物的百分比）为21.8%。各位点等位基因数介于2-9之间，各等位基因的大小范围在111-334bp之间；各位点期望杂合度（He）在0.1905-1.0000之间，各位点的观测杂合度（Ho）在0.1765-0.8090之间，除个别位点外（X781、X63、Z339）期望杂合度（He）与观测杂合度（Ho）均在0.3以上；各等位基因的多态信息含量（PIC）在0.1574-0.8272之间，其中X781与Z339两个位点的PIC值在0.25以下，属于低度多态性位点；X195等5个位点的PIC值在0.25与0.5之间，属于中度多态性位点；其余的17个位点的PIC值皆大于0.5，属于高度多态性位点。经Bonferroni法校正（P<0.003）后，24个多态性微卫星位点中有9个微卫星位点偏离了Hardy–Weinberg平衡，其中有2个位点（X1214、Z167）表现出显著的杂合子缺失（Hardy–Weinberg平衡偏离指数D<0，D=(Ho-He)/He），其他7个位点表现出显著的杂合子过剩（Hardy–Weinberg平衡偏离指数D>0）。（2）ISSR-PCR法开发微卫星标记。设计的8对简并引物中用于扩增CA/GT重复单元的引物PCA6、PGT6、PTG6、PAC6能够扩增出亮度较大的弥散带，将这四对简并引物克隆后测序。将前人已经开发出的微卫星同源序列去除后，共获得含微卫星的序列114个，成功设计引物45对。多态性分析显示，23对引物扩增条带为多态，开发效率为51.1%。各位点等位基因数介于2-12之间；各等位基因的大小范围在141-324bp之间；各位点期望杂合度（He）在0.2357-1.0000之间，各位点的观测杂合度（Ho）在0.3203-0.9169之间；各等位基因的多态信息含量（PIC）在0.2894-0.8877之间，其中4个位点（B78、B96、B75、D30）PIC值在0.25-0.5之间，属于中度多态位点，其他位点的PIC值都大于0.5，属于高度多态位点。经Bonferroni法校正(P<0.003)后，23个多态性位点中6个微卫星位点偏离了Hardy–Weinberg平衡，其中有三个位点（B96、B73、B92）表现出显著的杂合子缺失（Hardy–Weinberg平衡偏离指数D<0，D=(Ho-He)/He），另外三个位点（A89、B75、C2）表现出显著的杂合子过剩（Hardy–Weinberg平衡偏离指数D>0）。（3）微卫星序列特征分析对391条含微卫星重复的序列进行统计分析（图3-7），发现获得的微卫星序列长度范围在60-852bp之间，100bp以下的最少（16个）。大部分克隆片段大小在200bp以上，占总序列的83.1%。日本沼虾的微卫星序列中以二核苷酸重复最为丰富，占所有微卫星序列类型数目的92.4%，其次是三核苷酸（1.8%）或四核苷酸重复(1.6%)。二核苷酸重复中以（CA/GT）n类型所占比例最高，占总微卫星序列数的68.3%，其次是（GA/CT）n型，占25.1%，此外还有少量的（TA/AT）n型和（GC/CG）n型，分别占2.8%和0.6%。2.4个日本沼虾群体遗传结构分析利用8个微卫星标记（A28、A73、B13、B29、B44、B93、D7、Z337）对太湖、东平湖、微山湖、长江下游（崇明岛附近水域）4个群体的遗传结构进行分析，4个群体中长江下游（崇明岛附近水域）群体的平均等位基因数最多，为13.5000，其次是微山湖群体为12.2500，太湖群体最低为11.87。微山湖群体的平均期望杂合度和观测杂合度最高，分别为0.8403和0.8350；长江下游（崇明岛附近水域）群体的最低，分别为0.7145和0.7909。各群体间遗传距离分析显示，微山湖和东平湖群体的遗传距离最小，为0.0158；长江下游（崇明岛附近水域）和微山湖的遗传距离最大为0.3844。各群体间的遗传相似度分析也得出了类似的结果，微山湖和东平湖的遗传相似度最大（0.9162）；长江下游（崇明岛附近水域）和微山湖的遗传相似度最小（0.6312）。基于遗传距离的NJ和UPGMA聚类分析得到了相同的结果，东平湖群体和微山湖群体遗传距离最小首先聚类为一支，再与太湖群体聚为一支，最后和长江下游（崇明岛附近水域）群体聚类为一支。聚类分析的情况和实际的地理距离相符。4个日本沼虾群体的分子方差分析（AMOVA）显示，群体中89.29%的变异来源于个体内，仅有4.18%的变异来源于群体间，6.52%来源于群体内个体间，显著性检验分析显示这4个群体遗传分化极显著（P<0.01）。4个群体两两间遗传分化指数Fst从0.0031到0.0739，其中东平湖群体与长江下游（崇明岛附近水域）群体、微山湖群体与长江下游（崇明岛附近水域）群体的0.05<Fst<0.15，属于中等分化程度；其他群体间的Fst值均<0.05，分化程度较弱。日本沼虾群体的遗传结构分析对其种质资源保护和合理挖掘利用具有重要意义。（3）性别相关标记的筛选：在利用15雌和15雄日本沼虾对开发的47个微卫星标记进行多态性分析时，我们发现有8个SSR标记出现了雌性或雄性特异的等位基因。为此，我们合成了这8个微卫星标记的上游荧光引物，利用荧光SSR法，对另外的15雌和15雄日本沼虾进一步进行分析、验证。结果显示，有7个微卫星标记检测到了雌性或雄性特异的等位基因（等位基因出现次数在2次以上）。首先利用2个雌性个体和2个雄性个体对64对AFLP引物的扩增情况进行了筛选，对其中的9对引物进一步利用15个雌性个体和15个雄性个体进行大量分析。有8对检测到分布频率雌雄差异显著的等位基因，其中雄性占绝大多数（90.9%）的等位基因有1个，雌性占绝大多数（84.6%）的等位基因有2个，可为日本沼虾性别决定机制研究提供参考。
【作者】赵燕；
【导师】曾勇庆；王慧；
【作者基本信息】山东农业大学，动物遗传育种与繁殖，2014，博士
【关键词】日本沼虾；微卫星；AFLP；遗传结构；性别相关分子标记；

【参考文献】
[1]白如江,王尊新,鲍翠梅.智能化企业竞争情报系统研究[J].情报杂志,2004,08:84-85.
[2]刘一凡.我国环保企业发展现状与对策研究[D].北京交通大学,工商管理（专业学位）,2013,硕士.
[3]肖亚丽.国营林场森林资源信息管理专用软件开发研究[D].河南农业大学,森林经理学,2004,硕士.
[4]严梁恒.9种鸭科鸟类的线粒体基因组序列分析与雁形目系统发生研究[D].安徽大学,生态学,2014,硕士.
[5]赵婷.现代城市公共空间人工水景绿色设计策略研究[D].湖南师范大学,设计艺术学,2013,硕士.
[6]刘岱.美国报业体育新闻黄金时代（1920年代）历史与写作手法分析[D].北京体育大学,体育人文社会学,2013,硕士.
[7]黄布毅.DNA计算中若干理论问题的研究[D].华中科技大学,2005.
[8]傅强.石墨烯限域催化反应的原位研究（英文）[A].中国化学会.中国化学会第29届学术年会摘要集——第32分会：纳米表征与检测技术[C].中国化学会:,2014:1.
[9]刘雪华,邵鑫,刘凤珍,孟宪华,王敏,钟慧玲,宋金魁.氧化石墨-氧化铜复合材料的制备及其光催化性能研究[J].聊城大学学报(自然科学版),2013,01:40-42.
[10]姜敏.脐带缠绕的临床分析[D].吉林大学,临床医学,2014,硕士.
[11]龚涛.大直径盾构穿越地铁车站暗挖风道三维有限元分析[D].北京建筑大学,建筑与土木工程,2013,硕士.
[12]刘岚.WebNewsHunter智能代理[D].中国科学院研究生院（软件研究所）,2003.
[13]陈海.木结构建筑施工安全因素识别与安全行为模拟[D].厦门大学,建筑与土木工程,2014,硕士.
[14]谢诗琳.农村中小学教师聘任制改革研究[D].江西师范大学,教育管理,2003,硕士.
[15]邱陈辉.基于FPGA的自适应FIR滤波器设计与优化[D].江苏科技大学,信号与信息处理,2014,硕士.
[16]石丽琴.论信息化与可持续发展[D].广西大学,科学技术哲学,2004,硕士.
[17]樊世英1,2,3,孙军勇1,2,3,谢广发3,4,陆健1,2,3,5*.澄清剂对黄酒混浊蛋白去除效果的研究[J].食品工业科技.
[18]胡杰武,张秋生,胡靓.区域产业整合及经济增长研究-基于跨区域并购视角[J].中国软科学,2012,06:167-174.
[19]张伟.重结晶法制备细粒度HNIW[D].中北大学,安全技术及工程,2014,硕士.
[20]钟隽迪.回族伊斯兰教音乐研究中文文献综述[J].歌海,2014,02:19-23.
[21]陈莹秀.卡特福德翻译转换理论在农业科技英语汉译中的应用[D].南京农业大学,英语笔译,2013,硕士.
[22]曲阜阳.DDR SDRAM接口设计与静态时序分析[D].西安电子科技大学,集成电路工程,2012,硕士.
[23]寿文池.BIM环境下的工程项目管理协同机制研究[D].重庆大学,管理科学与工程,2014,硕士.
[24]罗人述.列车临时限速信息管理系统的设计与实现[D].广东工业大学,模式识别与智能系统,2014,硕士.
[25]董连平.具有水跃复锥结构的煤泥旋流重选柱研究[D].太原理工大学,矿物加工工程,2014,博士.
[26]刘莉.益髓生血颗粒治疗α-HbH病对造血细胞因子活性与表达的影响[D].北京中医药大学,中西医结合基础,2013,硕士.
[27]张倩.基于FEAF的离散型企业信息系统集成研究[D].成都理工大学,工业工程,2013,硕士.
[28]赵坤.内部控制质量与企业价值相关性的实证研究[D].沈阳理工大学,会计学,2012,硕士.
[29]孟婷.舞蹈啦啦队转体类难度训练方法的研究[D].华东师范大学,体育教育训练学,2013,硕士.
[30]刘海霞.X市S区城投集团公司战略转型研究[D].内蒙古大学,工商管理,2013,硕士.
[31]刘晓光.网络化制造中Web服务自动组合的若干关键技术研究[D].上海交通大学,2008.
[32]叶茂伟.封闭条件下风载荷对大型发射塔架动力学特性影响分析[D].哈尔滨工业大学,飞行器设计,2013,硕士.
[33]季淑娟,王瑜,张迎芳.数字化信息资源的整合与信息挖掘[J].冶金信息导刊,2004,03:32-34.
[34]刘璇.健身秧歌对老年肥胖女性身体素质的影响[D].山东体育学院,体育教学,2012,硕士.
[35]郑元润,张新时.毛乌素沙地高效生态经济复合系统诊断与优化设计[J].植物生态学报,1998,03:71-77.
[36]陈天佑.统筹城乡发展视角下的大理州“农转城”问题研究[D].大理学院,马克思主义基本原理,2013,硕士.
[37]周新栋.中文文本分类的文档索引机制及分类模型的研究[D].国防科学技术大学,2004.
[38]冯彩云.近自然园林的研究及其植物群落评价指标体系的构建[D].中国林业科学研究院,园林植物与观赏园艺,2014,博士.
[39]于宁,何通能,王泽锴.基于DSP的三相逆变并网同步控制算法及实现[J].机电工程,2013,07:878-881.
[40]张素婉.英汉有稿同传与无稿同传中的顺句驱动策略研究[D].南昌大学,翻译（专业学位）,2014,硕士.
[41]戎益升.试论城乡规划信息服务体系的完善[D].苏州大学,公共管理（专业学位）,2013,硕士.
[42]黄筱.编组站单双向布置方案选型方法研究[D].北京交通大学,交通运输规划与管理,2013,硕士.
[43]邓雅晴.锰铈复合型非贵金属氧化物催化剂的制备、表征及CO氧化反应机理的研究[D].华东理工大学,化学工程,2014,硕士.
[44]王宇良.LDPE釜式反应器混合特性研究[D].浙江大学,2014.
[45]孙洋.新闻客户端完善用户体验之路径研究[D].河北大学,新闻学,2014,硕士.
[46]程金路,张承进.自适应有限拍感应电机电流控制器设计[J].控制与决策,2005,10:1097-1101+1106.
[47]胡中泽.黄瓜花叶病毒在豇豆上诱导HR的无毒基因及功能域研究[D].南京农业大学,植物病理学,2013,硕士.
[48]于兆卿,艾淑萍,周晓平.大型热模锻压力机方案设计[J].一重技术.1999(02)
[49]张丽祥.论购房政策调整与情势变更的适用[D].兰州大学,民商法（专业学位）,2013,硕士.
[50]常文辉.妊娠早期、中期胎盘及引产胎儿HBV感染状况的研究[D].第四军医大学,流行病与卫生统计学,2004,硕士.